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We study the effects of a superconducting condensate on holographic Fermi surfaces. With a 
suitable coupling between the fermion and the condensate, there are stable quasiparticles with a 
gap. We find some similarities with the phenomenology of the cuprates: in systems whose normal 
state is a non-Fermi liquid with no stable quasiparticles, a stable quasiparticle peak appears in the 
condensed phase. 
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I. INTRODUCTION 



The problem of what happens when a large number of 
interacting fermions get together remains interesting de- 
spite many decades of work. The sign problem obstructs 
a numerical solution, leaving us to do experiments or the- 
orize. The metallic states of such systems that are well- 
understood theoretically are Fermi liquids. The basic 
assumption of this theory is that the states of the inter- 
acting system can be usefully put in correspondence with 
those of a collection of the same number of free fermions; 
in particular this means that the low-lying excitations of 
the system arc long-lived quasiparticles. 

This assumption fails in many strongly-correlated ma- 
terials. Quite a bit of effort has been made to understand 
what replaces the Fermi liquid theoryin the absence of 
stable quasiparticles [H, i, i i, i, 10, i, i, M M M 



fisl . ■ We believe that it is fair to say that it would be 
valuable to have a non-perturbative description of such a 
state of matter, hispired by work of Sung-Sik Le e llSl . a 
class of non- Fermi liquids was recently found [H, Il7l| (see 
also [l^, [l^ ) using holographic duality . This allows 
us to study observables of the strongly-coupled system 
using simple gravity calculations. For a review of these 

techniques in the present context, see 0, HI, ail. 

The analysis of applied to CFTs with a gravity 

dual, a conserved U{1) current, and a charged fermionic 
operator. Depending on the charge and dimension of the 
operator, it is possible to find Fermi liquid behavior, in 
the sense that the spectral function exhibits stable quasi- 
particles, or non-Fermi liquid behavior. At the boundary 
between these behaviors, one finds a marginal Fermi liq- 
uid, which arises as a phenomenological model [25j of the 
strange metal phase of the cuprate superconductors (the 
resistive state at temperatures larger than the critical 
temperature Tc for superconductivity, at a doping level 
which maximizes Tc). In this case, the contribution of 
such a Fermi surface to the resistivity also has the linear 
temperature dependence observed in the strange metal 

The calculation of the fermion spectral functions was 
done by solving the Dirac equation in a charged black 
hole background. The extremal Anti-de Sitter (AdS) 
Reissner-Nordstrom black hole (hereafter referred to as 
'RN'), which re pre sents the groundstate of the simple sys- 
tem studied in [l6l. Il7|. has a 'residual' zero-temperature 
entropy. This degeneracy is exact in the classical N oo 
limit; at finite N one expects it to be lifted to a large low- 
lying density of states. It is likely that the non Fermi 
liquid behavior does not depend on the large low-energy 
density of states: the small-frequency behavior depended 
on the existence of the IR CFT, not on the large central 
charge c oc s(T = 0) of the IR CFT. 

A closely related question regards the stability of the 
extremal black hole geometry. It is a stable solution 
of the Einstein-Maxwell theory. However, many known 
AdS string vacua which UV-complete this model con- 
tain charged boson fields which at finite density and low 
temperature will exhibit the holographic superconductor 
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instability [23, ■ Conveniently, the physical systems to 
which we would like to apply these models also generi- 
cally exhibit a superconducting instability {e.g. [29l.[30l|): 
the T ~ limit of most known non-Fermi liquids is un- 
der a superconducting dome^. The calculations in the 
RN black hole provide a model for the "normal" state 
above the superconducting Tc. 

As discussed in the last section of [13], this raises the 
following very natural question: what happens to the 
holographic Fermi surface in the presence of supercon- 
ductivity? One might expect to see a gap in the spectral 
weight, and we will see below that this is realized. Unlike 
the fermion two-point function calculation, here there are 
some choices for the bulk action. In addition to choosing 
the self-couplings of the bulk scalar ip, one must decide 
how to couple the scalars to the bulk spinor field C- It is 
always possible to include a |</'PCC couphng. In duals of 
matrix-like theories, where the spinor field is dual to an 
operator of the form tr A, it is natural to include a scalar 
(f with twice the charge of C, dual to the operator tr AA 
[3J]. Its dimension at strong coupling is not determined 
by this information. In this case, a (as it turns out, much 
more interesting) coupling of the schematic form CCv* is 
permitted by gauge invariance. We will specify the spinor 
structure of the coupling below. 

The effect of this coupling is to pair up modes at 
the Fermi surface, in a manner extremely similar to 
the Bogoliubov-deGennes understanding of charge exci- 
tations of a BCS superconductor. 

Interestingly, if the mass-to-charge ratio of bulk scalars 
is large enough, they do not condense [H] , and we pause 
here to comment on this case. This in itself is an inter- 
esting phenomenon which does not happen at weak cou- 
pling, and should be explored further. It means that the 
criteria for a string vacuum which exhibits the Fermi sur- 
faces described in [IB, [13 , but not the superconducting 
instability, are reminiscent of those required of a string 
vacuum which could be that of our universe: one doesn't 
want light scalar fields. In the latter context, a large ma- 
chinery [sl] has been developed to meet the stated goal, 
and similar techni ques will be useful here. In such a case, 
the calculation of [l6l . [l3| is valid to very low tempera- 
tures. One effect which cuts this off is the following^. In 
the RN black hole background, there is a finite density of 
fermions in the bulk [26| . There is a Fermi surface (in the 
sense that the bulk-to-bulk fermion spectral density has a 
nonanalyticity at uj — 0^ k — kp). There are interactions 
between these bulk fermions mediated by fluctuations of 
the metric and gauge field. The Coulomb force is naively 
always stronger [33|, but can be screened. This leaves 
behind the interactions by gravity, which are universally 



^ Other possible groundstates for holographic finite-density sys- 
tems, for example resulting from the presence of neutral bulk 
scalars, have been explored recently in |31|.|32|. Issll . 

^ We thank Nabil Iqbal for an instructive conversation on this 
point. 



attractive. There is some similarity with phonons. Of 
course, these interactions are suppressed by powers of 
N"^ (where TV^ = G^^ in units of the AdS radius). This 
may lead to BCS pairing with an energy scale 

T, ~ eY^e-^ - /ie-^' (1) 

where z^(0) is the density of states at the bulk Fermi sur- 
face, and V ^ N~'^ is the strength of the attractive inter- 
action. This is a very small temperature. This is exactly 
the scale of the splitting between the degenerate ground- 
states over which the RN black hole averages which is to 
be expected at finite N . Nevertheless, this is one way in 
which the RN black hole groundstate of the system stud- 
ied in [l^, [l3| is unstable, without the addition of extra 
scalar degrees of freedom. 

In this paper, we will probe (a few examples of) holo- 
graphic superconducting groundstates with fermionic op- 
erators. The retarded Green's functions Gii{uj,k) we 
compute may be compared with data from angle-resolved 
photoemission experiments on cuprate superconductors 
[SSl . [3^. In these experiments, a high-energy photon 
knocks an electron out of the sample, which is then de- 
tected. Knowing the energy and momentum of the inci- 
dent photon and measuring the energy and momentum of 
the detected electron allows one to infer that the sample 
has an electronic excitation specified by their difference; 
the intensity of the signal is proportional to the density 
of such states, A{LLi,k) = Im Gi? (cj, fc) (at least in the 
sudden approximation, which is believed to be valid for 
the relevant photon frequencies). Actual photoemission 
experiments have the limitation that they can only kick 
electrons out of occupied states, and hence can only mea- 
sure an intensity / oc A{u! ^ k) f (lo) , where /(w) is the 
Fermi factor, which at zero temperature vanishes for uj 
above the chemical potential. We do not have this limi- 
tation. 

Lest the reader get the wrong idea, we emphasize here 
several features of our calculations which differ from the 
experimental situation in any strongly-correlated elec- 
tron system. Perhaps most glaringly, as in previous work, 
the Fermi surfaces we discuss in this paper arc round. 
There is no lattice in our system. At short distances, 
our theory approaches a relativistic conformal field the- 
ory; the UV conformal symmetry is broken explicitly by 
finite chemical potential (we will also comment on the ef- 
fects of a small temperature). Also, our superconducting 
order parameter has s-wave symmetry, and so there are 
no nodes at which the gap goes to zero. It would be very 
interesting to improve upon this situation. 

Above the superconducting critical temperature T^, 
one usually has gapless excitations at k = kp. When one 
cools the superconductor below the locus {k = kp} 
generally remains the surface of minimum gap, i.e. the 
locus in momentum space where the gap in the fermion 
spectral density is smallest. This is not precisely the case 
here. This is because in general the holographic super- 
conducting condensate also affects the geometry outside 
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the horizon region, i.e. UV physics, and changes the ef- 
fective Schrodinger potential determining value of k at 
which the Dirac bound state occurs. The difference be- 
tween kp without the condensate and the surface of min- 
imum gap will be small in the examples we study, which 
have Tc small compared to ^, and are therefore close to 
the RN geometry as we review below. 

[ioj appeared as this paper was being completed. The 
crucial Majorana coupling is not included there. Related 
work will appear in [4ll . |42| | . 



A. Why the Majorana couphng is important 

We focus here on the case of odd d (the number 
of spacetime dimensions of the boundary field theory), 
where a single Dirac spinor in the bulk describes a sin- 
gle Dirac spinor operator in the boundary theory. In the 
case of even d, we will need to couple together two bulk 
Dirac fields. 

The bulk action we consider for the fermion is 

S[C] - J d'^+'xV^ [zC {T^'Dm - mc) C 

+ r^lv*eCT^(: + r^,v~^CT^^] . (2) 

tp is the scalar field whose condensation spontaneous 
breaks the U{1) symmetry. C is the charge conjugation 
matrix, which we specify below, and is the chirality 
matrix, {r^,r*^} = 0. The derivative D contains the 
coupling to both the spin connection and the gauge field 
Dm = Om + \i^MAB^^^ — iqc^M- We will occasionally 
refer to the coupling to the scalar in ([2]) as a 'Majorana 
coupling' because ('^CT^( is like a Majorana mass term. 
One reason for the necessity of the antisymmetric charge 
conjugation operator between the fermion fields in this 
term is that the simpler- looking object CaCa is zero be- 
cause the components are grassmann-valued. 

As we will describe, the coupling (p*(^'^C(^ + h.c. is also 
possible, but does not accomplish the desired effect. The 
coupling with the arises in descriptions of fermionic 
excitations of color superconductors |43j . In that context, 
the chirality matrix is required by parity conservation; 
since (p there is a bilinear of the same quarks to which it 
is coupling, the intrinsic parity of the quarks cancels out. 

One could worry that the perturbations of the scalar 
field will mix (in the sense that one will source the other) 
with the fermion equations of motion. This does not hap- 
pen in the computation of two-point functions because of 
fermion number conservation. 

We pause here to note the instructive similarity be- 
tween ([2]) and the action governing electrons in a BCS 
s-wave superconductor 



where a =T,i are spin indices, = I'^dfel — kp), and 
uj is measured from the chemical potential. This similar- 
ity is instructive because it explains why other couplings 
between the spinor and the condensate do not automat- 
ically produce a gap. 

The basis of modes which diagonalizes such an action 
is the Nambu-Gork'ov basis: 



7„(fc) = u{k)c^{k) + C1v{k)c%-k) 



(4) 



note that u and v do not have spin indices. The Green's 
function which results from this mixing is 



(cfe(w)^Cfc(tj))j 



■a 



(w -I- ie) 



f2 ■ 



|A(fc)P 



(5) 



This function has two poles for each fc; they approach 
Re {lo) = as fe approaches the Fermi surface. Each has a 
minimum real part of order A. The residues of these two 
poles, however, varies with k: at large negative k — kp, 
the weight is mostly in the pole with Re (w) < and the 
excitations is mostly a hole. As k moves through kp , the 
weight is transferred to the other pole, and the excitation 
is mostly an electron. Without such a mixing between 
positive and negative frequencies, the Green's function 
would have only one pole, which would be forced to cross 
Re {uj) = as /c goes from k <^ kp to k ^ kp, and there 
could not be a gap. This continuity argument assumes 
that in the absence of the condensate, the dispersion is 
monotonic. 



II. REVIEW OF GROUNDSTATES OF 
HOLOGRAPHIC SUPERCONDUCTORS 



Consider the action 



^=^[R+J2- lidA)' - |(V - ^q^A)^\' - ml\^{- 



(6) 

We will work in units where the AdS radius L is unity. 
For — 2q^ < —3/2, the Reissner-Nordstrom AdS solu- 
tion is unstable at low temperature to forming scalar hair. 
The extremal limit of these hairy black holes was found in 
[4^'^. Unlike the extreme Reissner-Nordstrom black hole, 
the area of the horizon goes to zero in this limit. The de- 
tailed behavior near the horizon depends on and q^, 
but for TO^ < 0, the solution has Poincare symmetry 
near the horizon. This has an important consequence. 
Consider solutions of the Dirac equation with e*'^''^ de- 
pendence. If k is timelike in the near horizon region, then 
one can impose the usual ingoing wave boundary condi- 
tion to compute the retarded Green's function Gr. Since 



S[c] = j d^-^kduj {ci{uj, k) {iio - 4) c^{lo, k) (3) 
-A(fc)c|(tJ, fc)c|(-w, -k) - A*(fc)c|(w, fc)cx(-w, -k) 



^ Groundstates of holographic superconductors, inclu ding other 
forms of the scalar potential, were also studied in |45|I . Our 
analysis should apply to those whose IR region is AdS a \ we leave 
the other cases for future work. 
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the boundary condition is complex, the Green's function 
is complex, and hence Im Gjj is nonzero indicating a 
continuous spectrum of states. However, if k is spacelike, 
the solutions are exponentially growing or damped. Nor- 
malizablility requires the exponentially damped solution. 
This is a real boundary condition, and so the solutions 
will be real and Im Gr = 0. This is qualitatively dif- 
ferent from the extreme Reissner-Nordstrom AdS whose 
near horizon geometry is AdS2 x R^. In that case, there 
is a continuous spectrum for all (w, /c'). 

The light cone in the near horizon region will not have 
the same speed of light as the asymptotic geometry. One 
can show that as one approaches m^ — 2gJ = —3/2 where 
the RN solution becomes stable, the speed of light in the 
IR CFT, ciR (not to be confused with the central charge 
of the infrared CFT), goes to zero (sec FIG. [T|). This 
means that in momentum space, the light cone opens up 
so all momenta are effectively timelike, and the spectrum 
continuously matches onto the RNAdS case. 
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2(2a2 + 7a + 5) 



+ 5q: + 6 ~ 



1/2 



(11) 

Although the curvature remains finite, derivatives of the 
curvature diverge at r = unless a = 0. FIG. 2 shows 
the solution for g(r) and (j){r) for a choice of q^p which is 
close to the value •\/3/2 where Reissner-Nordstrom AdS 
is stable. One sees that g dips down and has a local 
minimum at a value r « 1. As ^3/2, g vanishes 

at this local minimum which becomes the horizon of the 
extremal Reissner-Nordstrom AdS black hole. 

For TO^ < ( and > — m^/6), the zero temperature 
solution near the horizon is 



= 2(-logr)i/2, g^{2mlliy\ogr, 



= 0or^(-logr)i/2 



where 




-K\ogr 
(12) 

(13) 



(14) 



and 4>o is adjusted to satisfy the boundary condition at 
infinity. The near horizon metric is (after rescaling t) 



2m^r2 log? 



(15) 



FIG. 1: The speed of light in the IR CFT, cm, as a function 
of the boson charge. The blue thick curve is = — 1, the 



red thin curve is m„ 



0. The vertical dashed lines indicate 



the value of below which the RN solution is stable. 



In more detail, the static, plane symmetric solutions 
take the form: 

dr^ 

ds^ = -g{r)e-^^'-Ue + — - + r''{dx' + dy'') (7) 



One clearly sees the Poincare symmetry (but not the con- 
formal symmetry) in this case. There is a rather mild null 
curvature singularity at r = 0. 
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A = (j){r) dt, If = Lpir) 



(8) 



For mfp = 0, the zero temperature solution not only has 
Poincare symmetry but approaches AdS^ near the hori- 
zon, and r = is just a Poincare horizon. The leading 
order corrections can be found analytically and depend 
on a parameter a which is a function of q^p, but stays 
small (|a| < .3). Explicitly, 



(9) 



FIG. 2: This plot of the emblackening factor g (left) and 
the electrostatic potential (j) (right) in the — 1.3, ni^ — 
groundstate solution exhibits the almost-RN horizon at r = 1. 
In this plot and those below, we use units where fi = 2\/3. 



III. DIRAC EQUATION 

The Dirac action is 



where 



0^ + 5(1 + 6 



1/2 



Xi 



+ 5a + 6 
4(a+l) ' 



(10) 



So^i I d'^+'x C {r'^DM - - A|</.|2) C (16) 



where we are using the conventions of [46|. The A cou- 
pling could be replaced by a more general function of 
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We will set A = for now except to discuss its 
effects briefly below. 

As discussed in section II Al if the charge of the scalar 
is such that = 2qQ then we can add to this 



Sr. = / d''+ 



-5^*e(^* + %*r5)c + h.c (17) 



where the charge conjugation matrix is 

= cr^c* (cr^) (cr^-) 



(18) 



This term is essentially a Majorana mass term. There are 
two terms because there are two Majorana spinors in the 
bulk (or Weyl spinors) and these can have independent 
masses. 

In the case of odd numbers of bulk dimensions, there 
is no and this term does not exist. This matches the 
fact that in odd numbers of bulk dimensions, a single 
Dirac spinor in the bulk describes a chiral fcrmion oper- 
ator in the boundary theory (47j : such a fcrmion cannot 
be paired with itself in a rotation-invariant way. The 
analogous coupling in odd bulk dimensions requires two 
Dirac fermions. That this is possible can be seen by di- 
mensionally reducing a theory with an even number of 
bulk dimensions on a circle. We will not discuss this in 
detail here. 

Now we study the Dirac equation in more detail. It 
turns out that the same simplification that appeared in 
the RN background occurs for the more general metric 
([7]). Very briefly, the form of the spin connection 



uJif = dt ^/g^^r {yfgtt) wj, 
implies that 



(19) 



ic.,,Mef r^r'^" = h^dr In (-gg-) (20) 

so we can rescale J- = {—gg^'^Y^^C ^-nd remove the spin 
connection completely. The new action is 



So 



i J d'^+^Xy^T {T^'D'j^j - mc) T 



(21) 



where D'j^j ~ dm ~ ig(A]^J with no appearance of the spin 
connection. 

The Dirac equation following from 5*0 + S",; is 

-mc_)T + 2iip{T] - r]5T'^)CT*-T* = . (22) 

Expand this into Fourier modes with = k, ky = 0: 

{Ip'{k, uj) - mQ)T{k, Lu)+2iip{rj-r]5r^)CT^T*{-k, -lu) = 

(23) 

To get any further we must specify a basis of Dirac ma- 
trices. We focus on d = 3, that is, a 3 -I- 1 dimensional 
bulk. We choose a basis of bulk Gamma matrices as in 
0, 



-cr3 
-a' 



-a^ 
a2 



pa 









ia' 
-ia^ 



(24) 



such that F-* = — F- and F-* = F- which fixes the charge 
conjugation matrix to be CF- = F-. This basis has the 
features that (with 775 = and ky = 0) the Dirac equation 
is completely real. 

We will now split the 4-component spinors into two 2- 
component spinors = {^1,^2)'^ where the index a = 
1, 2 is the Dirac index of the boundary theory. Then 



= i^Dr{±k) + Vr*cr^'^j J^i,2(fc,w) 

~2ic7^ipriTl2{~k, -uj) ± 2icrV%^2* (25) 
where^ 

Dr{k) = -VF^cr^S,. - - ^f^ia'^k -\- y^a^q^At 

(26) 

We see that the 775 term mixes J^i{k, lu) with J^|(— /c, —lu) 
(and J-2{k, lu) with J-^{—k, —uj)) - this is the mixing that 
will most interest us, because for the RN background 
these two fields have coincident Fermi surfaces (at a; = 0). 
Setting 7/ = ([25|) becomes 



Dr{±k) 
where 



g^*-LU 



±2i(pri5 



Ti{k,uj) 
J^i{-k,-Lu] 



±2iip*rj; 



*2 



1'l,2 = 



J^2ik,Lu) 

Tl{-k,-u] 



(27) 
(28) 



are the bulk analogs of the Nambu-Gork'ov spinor.'^ We 
see explicitly from ((27|) that for a general black hole back- 
ground in the absence of mixing (775 = 0) the spectrum 
of Ti{k,Lu) compared to the spectrum of J-"2(— fc, —lu) is 
a refiection about the cu = axis. This is crucial for 
generically generating gapped states for non-zero 7/5. 

We have set 77 = both to make the analysis easier 
and because turning on both rj and 775 implies that some 
discrete symmetry of the boundary theory is broken. 

For completeness, we record the Dirac equation with 
775 = 0, 77 7^ 0. In this case, the mixing would be between 
J-i{k,Lu) and J-l(—k, —lu) with the equation being 



Dri±k) 
DriTk) 



'g*-*-Lua^ —2iipria^ 
g**Lua^ 



2itp*7j*a-^ 



*1,2 = 

(29) 

where ^1^2 = {^i,2{k,Lu),J^i 2i^k,—Lu))'^. Because the 
differential operators Dj. in the diagonal entries above 



^ The frequency uj is measured from the chemical potentiaL 
^ The index on „ is the boundary theory Dirac index. For the 
rest of this section l llllll for simpHcity of the discussion we will 
concentrate mostly on one of these: N&i. In section ||IV[| wc will 
give results for ^2- 
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are evaluated with opposite fc, the two mixed compo- 
nents will not have coincident spectra at a; = 0, 77 = 
(see Figure 5 of to see this in the RN background). 
As such there will only be eigenvalue repulsion if there is 
some accidental eigenvalue crossing, and this will generi- 
cally occur away from uj ~ 0. This should be contrasted 
with the 775 mixing discussed above. 

A. Boundary conditions 

As reviewed in section 2, many of the solutions found 
in [4^ have an emergent Poincare symmetry in the deep 
IR, and some even have emergent conformal symmetry. 
For now we will mainly consider the latter case in which 
the solution approaches ^^54 near the horizon. To de- 
termine the IR boundary conditions for the spinor ap- 
propriate for the retarded Green's hmction, wc consider 
the Dirac equation in the far IR region, where the metric 
is just pure AdS^ with no electric field and zero chemical 
potential: 

ds^ = r2 {-cf^dt^ + dx^) + ^^i^ 
= (fi^Lpo x = Xo- (30) 

The speed of light in the dual IR CFT is cir = 
e"^"/^ / Lju. The most relevant terms in the Dirac equa- 
tion close to the Poincare horizon are 9,.'I'i = with 

Very generally, the off-diagonal terms are subdominant, 
by arguments given in [4J| in the discussion of the 
Schrodinger potential for the optical conductivity: the 
relative magnitude of the off-diagonal term to the terms 
appearing in (|5T|) is (p^/gu = ^\J1)&~^^'^ which must gen- 
erally vanish on the horizon. 

Because of the diagonal form of (|3T|) . wc can con- 
struct a basis of ingoing solutions by considering T\ (fc, u) 
and T^^—k^—Lo) separately. As is familiar from zero- 
temperature AdS^ the character of the boundary condi- 
tions depends on the sign of = —Lu^/cj^ + k'^. To begin 
with we will work outside the light cone where > is 
spacelike. Here the behavior of the solutions is normal- 
izablc and non-normalizable. Wc will pick the one which 
is normalizablc at r ^ 0: 

(I) T;\-k, -uj) "^S" 0, Tl{k,u:) 4e-^^-/^ = 

is an eigenvector of the matrix AI in (j3ip . This now 
allows us to formulate the general incoming boundary 
conditions in order to compute retarded correlators. Wc 



simply use the ie prescription to define how to continue 
the branch cuts in ((5^ to timelike < 0. That is, we 
take Lu ^ uj + ie. 

For the other component !F2{—k^—tjj) we can simply 
take ijj —ijj in (j32p . 

(II) Ti\k,ujy^\, .f*"(-fc,-c.)'^~°4v^^-A- = 

and again for timehke < in pS]) we continue lo — > 
uj + ie.^ In the absence of the 775 mixing, the two solutions 
of the Dirac equation (|27p for '^i specified by the IR 
behavior I, II compute the Green's functions for the two 
boundary fermion operators, 01(^1 k) and 0\{—uj, —k). 

Now wc consider the boundary conditions at the 
boundary of the UV ^^5*4. This will tell us how to read 
off the field theory correlators. The mixing term is again 
subdominant at the UV boundary, so that the asymptotic 
behavior is the same as usual: 

(34) 

and similarly for I II. The boundary retarded Green's 
function is: 

(B\ /G^.ol Go.o^l A\ All X 

\b*^ B*"; - l^G^tc^t G^tc^J \-Al^ -Afy 

(35) 

The definition of the Green's functions appearing above 
is: 

GcD{i^,k)^i j d'^-^xdte'''^~'"^'e{t)({C{x,t),D{0,0)}) 

(36) 

Note that the spectral densities (which should be positive 
by unitarity) arc ImG^tQ^ and ImG^^^^^t. 

More generally including the analysis for '^'2 the above 
matrix (j35p will fit into the Lorentz covariant correlator 
which is a 4 X 4 matrix (recall that this is for k^ ~ k,ky = 



° Beware the following confusion: because there is a complex con- 
jugation on 1 oris might expect this to switch the sign of the ie. 
This is not the case because we should think of analytically con- 
tinuing JF|(— fc, — a;) — > fc, — tj*); this procedure preserves 
the incoming boundary conditions. 

^ The minus signs appearing in front of come from the fact 

that — (Aj)^ is the source for where the minus sign is from 
anti-commuting this (Grassman valued) source in the boundary 
theory action so that it is in the correct order and the action is 
real. 
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(37) 

where O = (01,02)^ and Oc = (C7*)(Ct)'^ where the 
boundary theory charge conjugation matrix can be shown 
to be C7* = 1. Note that all the entries in this 4x4 
matrix will be non-zero if both 77, 775 are turned on. 



B. Evolution equation 

It turns out there is a super nice way to package the 
above linear differential equation into a non-linear evo- 
lution equation, in the spirit of the evolution equations 
considered in [43| and used in [1^, [Z^] . This is useful for 
identifying the Fermi surfaces numerically, because al- 
though the spinor components vary greatly with r, their 
ratios, which appear in the evolution equation, remain 
order one. 

Define the following matrices: 



Y 



Z = 



G = YZ- 



(38) 



Then one can write the following evolution equation: 



(39) 



G { ^ka^ + y^(c. + qcAa^) + ( 'J j ) G 

-V5 



V5 



+ -Vg^ka-' + V5"(w + Qc^tcr ) + 2itp 



The boundary conditions on this matrix at the IR 
AdSi horizon and the UV ^^54 boundary become re- 
spectively: 



G 



G 







(40) 



IV. RESULTS: BOUND STATES OUTSIDE THE 
EMERGENT LIGHT CONE 

A. No mixing 

We start by looking at 7/5 = = so there is no mixing. 
We will concentrate on the field J-2{k, ui) (from which we 
can reflect about w = to generate !Fl{—k, —lo). ) Note 
that we are now switching 1 2 relative to the discus- 
sion of the previous section - all results in this section 
will be for the Nambu Gork'ov spinor '^2- The reason 
being the primary Fermi surface in the RN background 
(the one with largest kp) makes its appearance in the 
Green's function for T2{k,uj) [l^- We are interested in 
understanding the fate of this primary Fermi surface in 
the condensed phase. 

Now since the initial conditions are real for spacelike 

> and the Dirac equation for T2 in the absence of 
mixing is real, the spectral functions should be zero out- 
side the emergent IR liglitcone. This is true up to delta 
functions which can appear because the real part of the 
Green's function has a pole which becomes a delta func- 
tion in the imaginary part thanks to Kramers-Kronig. 
These are bound states of the Dirac equation since they 
are normalizable at the IR AdS/i horizon and at the UV 
AdS/^ boundary. They will represent infinitely long lived 
fermion states in the field theory. 

For now we will look for these states in a small set 
of the zero temperature hairy black holes constructed in 
and reviewed above. We will concentrate on the case 
with zero scalar potential energy {V{ip) = ^ tt^^ = 0) 
with general charge for the scalar. In this case L/_r = 
1 and the speed of light in the IR CFT can be found 
numerically, sec FIG. [T] 

The fermion charge* will be constrained by gauge in- 
variance to be = q^p/'i, so that the 775 term can be 
added later. The mass of the fermion is a priori indepen- 
dent of the mass of the scalar. We work with mc^ = 
for numerical convenience. It will be interesting to look 
at small charges close to the critical charge q^p a/3/2 
where the critical temperature Tc ^ and cir, — > 0. 
For q^ < V3/2 the RN black hole is stable, and as can 
be seen from FIG. [21 the superconducting groundstate 
approaches the RN solution. In this limit the spectral 
densities should look more and more similar to the ones 
of the RN black hole, which we have a good handle on. 
Indeed, for reference, wc know that there is a Fermi sur- 
face in the RN black hole for 777^ = and = •\/3/4 
when kp w 0.75 with IR scaling exponent v w .18. 

FIG.[3]shows the location of these states for different q^ 
in these zero temperature superconducting backgrounds. 



Note that when 775 = the evolution equation preserves 
the diagonal form of the initial condition in the IR. 

This method runs into difficulty if Z becomes non- 
invertiblc at finite r; this happens for the multi-node 
boundstatcs associated with secondary Fermi surfaces. 



There is a faetor of two difference in the normalization of the 
charges for both scalars and spinors in [Tsl (LMV) compared to 
[3 (HR) - they are related by qlmv = "^QHR- Wo will work 
with the qufi normalization throughout. 
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FIG. 3: Boundstates outside the IR lightcone for various val- 
ues of q^p = 1.5, 2, 3, 5. Note that the frequency axis has been 
scaled by cir which is different for each curve. We can't resist 
mentioning the approximate relation kp « for where the 
curves cross uj = Q. 



B. Mixing 

The stable gapless (w = 0) excitations we have found 
in FIG. 3 seem rather surprising in a strongly coupled 
theory. We now demonstrate that turning on 775 7^ 
(and keeping ?7 = 0) the stable excitations studied above 
develop a gap. The reason for this can be simply under- 
stood by the general arguments of eigenvalue repulsion. 
Since the positive frequency modes mix with the nega- 
tive frequency modes (at the same k) the repulsion occurs 
when w = 0. 

More carefully, we can study the Dirac equation with 
mixing. Because the initial conditions ((5^ and for 
spacelike < are real one might expect that again 
the spectral functions are zero except for delta functions. 
This is a little subtle because the Dirac equation (P7)) is 
real except for the 775 term. However it turns out that 
despite this, the spectral functions are still zero. We 
can see this in two ways. Firstly, the spectral functions 
are the difference in the retarded and advanced Green's 
functions (this is more general than the imaginary part 
of the retarded function). For spacelike > these 
two Green's functions arc calculated with the same Dirac 
equation and the same initial conditions (the difference 
comes from the it prescription when going to < 0.) 
Hence Gr = Ga here and the spectral function is zero 
except for on bound states. 

Secondly, the evolution equation ([5^ for spacelike 
> preserves the following form of the 2x2 Green's 
function matrix G (recall we have switched 1 2 relative 
to ^ ): 

(41) 

Hence the spectral densities for G^^^^f , G^^^i arc zero. 
The phase of 7/5 is arbitrary since we can change it by 
rephasing the operator O, hence it cannot matter for the 
spectral density of GoiOi- 

To find the bound state in this situation we should 
look for places where det G~^ = at the boundary. Note 



that dctC""'^ G M for > so indeed this is a well de- 
fined problem. This delta function will appear in all 4 
spectral densities. The residue however will be different 
in each component. We concentrate on G^^^^t because 
this is what should be accessible to photoemission "ex- 
periments". The results are given in FIG. |3] and FIG.[5l 




0,0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

i k 



FIG. 4: Mixing between positive- and negative-frequency 
modes due to the Majorana coupling. Shown are density plots 
of the fermion spectral density A{k, cu) — ImG^ for g,; = 

|,mj = 0. The first plot is in the T = RN black hole, 
no scalar. The remaining plots are in the zero temperature 
background with — §,771^ — 0, for various values of the 
Majorana coupling, rj^ = 0,0.2, 1.5. 

We can learn something from perturbation theory in 
?75 . The splitting is determined by the eigenvalues of the 
matrix 

where 

P^^J drV^xL°^c.v7^xL"'(-l)" = ^JL (43) 
(J was defined in [l^ . appendix C) and 

Ql = y dry/g:;:;xf^2ir]5ipxf\Qi = J dr./g^Xi''2irilip*x'^ 

(44) 

where Xq*^ denotes the boundstate wavefunction in the 
basis X] = -^ij Xi = -^2 (""^i —k). Thinking of the Dirac 
equation as a Schrodinger problem, this matrix V is the 
perturbation Hamiltonian in the degenerate subspace. 

The fact that at w = 0, 775 = 0, the up and down 
boundstates are the same implies that P| = — P| = P 
and Q| = Q*); the eigenvalues of V arc therefore 

± V^P^+W- (45) 

Looking for low-energy boundstates with fixed k then 
requires these eigenvalues to vanish, which occurs when 
-P2 + |Q|2 = 0, i.e. when lo - \ri^\. 
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FIG. 6: The effective Schrodinger potentials for a probe 
scalar with mass mp^obe = —3/2 and gprobe = 5. The hori- 
zontal axis is the tortoise coordinate; the UV boundary is to 
the right, at f = 0. The different curves are different values of 
uj > 0; the top curve in each plot is for ui = 0. Left: The po- 
tentials for the groundstate with = and = 1. Right: 
The corresponding pictures for the RN black hole with the 
same charge density. 
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FIG. 5: The effect of the Majorana coupling on the fermion 
spectral density. Shown are plots of A{k,ui) at various k G 
[.81, .93] for qQ = ^,m(; = in a low-temperature background 
of a scalar with q^ = l,m^ = —1, with 7)5 = 0.025 (top) 
and rjs = 0.075 (bottom). The blue dashed line indicates the 
boundary of the region in which the incoherent part of the 
spectral density is completely suppressed, and the lifetime of 
the quasiparticle is infinite. The red dotted line indicates the 
location of the peak. 



C. Luttinger-like behavior near the lightcone 



To understand what's happening at u'^ 



we 



consider the Schrodinger form of the wave equation, 
where the role of the energy eigenvalue is played by — fc^. 
For simplicity (and because the pictures are nicer) we 
draw the potentials for the case of a charged scalar probe 
(not to be confused with the charged scalar ip which is 
condensing.) For further details, see Appendix B of 

The physics of the IR lightcone is visible in FIG. [D In 
the RN background (right plot), turning on any nonzero 
frequency opens up a bottomless pit in the effective po- 
tential leading into the AdS2 region where the tortoise 
coordinate f — > —00. Therefore, in the RN groundstate 
there are no infinitely-stable quasiparticles with nonzero 
frequency. On the other hand, in the superconducting 
groundstate, the limiting value of the effective potential 
as f — > —00 is —Lo'^/cjjf . Therefore, there is a threshold 



frequency = |ciRfc| below which the IR limit of the 
Schrodinger potential remains above the boundstate en- 
ergy. More precisely, there will be a normalizable bound 
state close to the boundary as long as the energy (— fc^) 
is less than the limiting value —oj'^/cjj^. Beyond this the 
bound state enters the light-cone and is no longer a stable 
quasi particle. 

The fact that we see a stable particle below the contin- 
uum is qualitatively what one expects for systems with a 
gap Wo- For energies ujq < lu < 2u!o, one excites a single 
quasiparticle which is stable since there is nothing for it 
to decay into. Only at energies above 2ujo does one start 
to see a continuum. 



The spectral density near the lightcone, and in particu- 
lar the width of the quasiparticle after it enters the light- 
cone can be computed by matching between the ^^5*4 
regions in UV and IR as in [13, SJ] . The size of the over- 
lap region is controlled by the quantity = k"^ — uP' j c^^i, 
which should be small in units of the chemical potential. 
In the notation of , the result for the Green's function 
is of the form 



(46) 



where A± , B± are reaf data associated with the UV re- 
gion, and Q is the IR CFT Green's function to be dis- 
cussed below. If there is mixing between positive and 
negative frequency modes then A± , B± are 2x2 ma- 
trices in the basis of the Nambu-Gork'ov spinor. They 
are smooth (analytic) functions of A;, w so the leading non 
analytic behavior in /c, w is from Q . For purposes of expo- 
sition we will describe the results for a probe scalar field 
in parallel to that of the spinor. We will leave details of 
the spinor calculation to Appendix A. 



® They are only real if we take r?5 G 
loss of generality. 



which wc can do without 
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For a probe scalar, the IR CFT Green's function is 







V 







(47) 



The quantities are related to the IR CFT scaling di- 
mension of the boundary operator by ^ffj = ^ + , and 
are determined by studying the behavior of the field at 
the UV boundary of the IR AdSi region in They 
are given by 




probe 



± l?75|vo) 



(48) 



where ip^ = ip{r = 0) (the subscript c is for 'condensed' 
and is intended to distinguish this object from the analo- 
gou s IR CFT scaling dimension in the AdS2 region of RN 
[17j). Notice that the IR CFT scaling dimension depends 
on the coupling rj^. 

For the probe spinor the IR CFT Green's function ap- 
pearing in (|46p is 








(49) 



For the spinor case the relation between Aj/j and 775 



IS, 



(50) 



see Appendix A for more details. 

We can extract two interesting statements from these 
calculations. From the form of (and in particular 
the reality of A, B) we learn that at generic w, k (but 
small |s| so that this matching applies). 



ImGcx (B_ - B+A-'^ A^) {liag) A:;'^ 



(51) 



The dependence of Vc on 7^5 has the following conse- 
quence. In the last plot of FIG. 21 one can see that 
the coupling to the condensate is also suppressing the 
incoherent spectral weight inside the lightcone. This is 
because the IR CFT dimension is becoming large as we 
make 775 large. 

Finally, if we look near a quasiparticle pole, which close 
to the light-cone occurs when det = 0, we see that the 
imaginary part of the location of the pole is determined 
by the IR CFT Green's function. This determines the 
width of the resonance as it enters the lightcone. The 
result is that the width behaves as 



(lU - CiRfc) 



(52) 



for the scalar and spinor respectively, which can be com- 
pared to the behavior in FIG.IH^'' 



Actually we need to be more careful for the case Uc < 1/2 (for the 
spinor.) Here 1521 1 should be replaced by, F ~ (oj — cjj^fc)^''^^'''^ . 



We emphasize that there are two mechanisms which 
suppress the spectral weight: one sets it exactly zero (ex- 
cept for delta functions) outside the light cone. This is a 
property of IR behavior of the background geometry. The 
other mechanism suppresses the weight independently of 
the momentum (this is a numerical observation visible 
from the dashed blue lines in FIG. 0, and depends on 
the scalar-spinor coupling. This mechanism generates the 
gap for the quasiparticle peak, and can be understood in 
terms of the dependence of the effective IR scaling di- 
mension of the fermion operator on 775 as in the previous 
discussion. The latter mechanism also affects physics at 
k = whereas the lightcone mechanism does not. 



V. DISCUSSION 

We should make a few remarks about the effects of 
other possible couplings between the bulk spinor and 
scalar. The coupling 



utral 



[C] 



-gXlvfCC (53) 



is possible whatever the charge of the spinor and scalar. 
By the argument given in section \I~K[ the Green's func- 
tion for the system with 775 = near kp should have only 
one pole (whose location may however be dramatically 
affected by the couplings A, 77), and the effects of the in- 
teraction (|53p cannot be interpreted as mixing of particle 
and hole states. As the |vPC^ coupling is varied, it is easy 
to be fooled into thinking that there is a gap even when 
there is not, when looking at energy distribution curves 
because the Fermi momentum moves with A. 

As observed first in ^l| increasing the mass of the 
effective field in the IR (which can be achieved by ei- 
ther including the above A coupling or changing the UV 
mass: mjn = + X(Pq) can lead to poles which never 
reach the w = axis and may be interpreted as gapped. 
This mechanism for removing low-energy spectral weight 
(which happens because increasing tti/a pushes up the 
effective Schrodinger potential for the Dirac equation) is 
qualitatively different from the mixing described in the 
previous sections. It is analogous to adding a relativistic 
mass to particles and anti-particles in a relativistic field 
theory at non-zero chemical potential. The gap in this 
situation is around uj ~ — /.( and does not generically pro- 
duce a gap at w = 0. This should be compared to the 
gap from the 775 coupling which is like adding a mass to 
particles and holes (absence of particles) about the Fermi 
surface at a; = 0. 

It was shown in (43 | that the zero temperature super- 
conductors we have studied here do not have a hard gap 
in the optical conductivity: The real part of the conduc- 
tivity remains nonzero (although typically exponentially 
small) at low frequency and T = 0. Despite this fact, 
one might have wondered whether such a hard gap in 
the conductivity exists for the fermionic probes that we 
study in this paper. The existence of a non-zero spectral 
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weight around the origin of FIG. |4] suggests that this is 
not the case; however, to see this effect in the conductiv- 
ity it would be necessary compute a X/N"^ correction as 
in d^. 

It would be interesting to understand better what 
property of the boundary theory is reflected by the pres- 
ence of the ?75 coupling, which is required to produce an 
actual gap in the fcrmion response. One clue is that its 
presence specifies the 'intrinsic parity' of the dual oper- 
ator, i.e. the dual operator acquires an interesting phase 
under a parity transformation. Realizing string vacua 
where this coupling is nonzero would probably be valu- 
able. 

So far we have considered the fermion spectral func- 
tion at zero temperature. FIG 7 shows what happens 
as one raises the temperature. The temperatures shown 
are much less than Tc. As T ^ Tc, the condensate goes 
to zero, so its coupling to the fcrmions goes to zero and 
the gap disappears. Actually, the thermal broadening of 
the peak makes the gap disappear at about .7Tc. In the 
opposite limit, as T — > 0, the width of the peak vanishes 
rapidly. It appears to vanish faster than a power law, but 
the general temperature dependence deserves further in- 
vestigation. 




0^ ^ Nr-^ I VV , -~^— ry^ i w^ ... I 

-0.010 -0.005 0.000 0.005 0.010 



FIG. 7: The effect of temperature (much less than Tc) on the 
fermion spectral function. Shown are plots at = l,m^ = 
—1,(7^ = |,m^ =0,775 = .025, and momenta where the peak 
is closest to tj = 0. The different curves correspond to differ- 
ent temperatures approaching T = Q. 

We close with a few comparisons with real phenomena. 
Here we make a simple observation which follows from 
the sharpness of the peaks in the 'no man's land' regime 
(i.e. outside the IR light cone). This regime is induced by 
the superconducting order. This means that if we start at 
high temperature in the normal phase with some Fermi 
surface without stable quasiparticles (like say a marginal 
Fermi liquid case, v ~ ^ 'm the notation of [13 )j and 
cool into the superconducting phase, sharp quasiparticlc 
peaks appear, at least for 775 not too big. This matches a 
mysterious piece of cupratc phenomenology: in the nor- 
mal phase, photoemission experiments show no stable 



quasiparticlc peak, but a coherent peak emerges in the 
superconducting phase (see e.g. figure 47 of the review 
[50|). From the gravity point of view, this is happening 
because the scalar condensate is removing the AdS2 re- 
gion which was responsible for the finite lifetime of the 
holographic quasiparticles [l3l: this is the gravity state- 
ment that the condensate is lifting the many gapless ex- 
citations into which the quasiparticlc could decay. The 
mechanism for the stability of these excitations is very 
similar to the recent holographic explanation [5l| of the 
critical velocity in a (holographic) superfluid below which 
there is no drag, and above which energy is dissipated by 
the creation of IR AdS^i unparticles. 

This similarity can be made more precise. In a BCS 
super/?Mid, the decay of the quasiparticles can be medi- 
ated by emission of a Goldstone boson (this mode is eaten 
in a superconductor, and the following effect is absent). 
It can happen that this decay is kinematically forbid- 
den: the decay cannot happen if the group velocity of 
the quasiparticlc is larger than the speed of sound (see 
appendix B of [s^). In our system, the quasiparticles 
develop a finite lifetime when they can decay into the 
modes of the IR CFT dual to the IR AdS^ region. These 
modes are distinct from the Goldstone mode (which is 
apparently hidden by powers of N), but the effect is the 
same. 



The energy distribution curves {A{k,uj) at fixed k) 
shown in FIG. O exhibit another feature in common 
with ARPES measurements on the cuprates, namely the 
so-called 'peak-dip-hump' structure: in addition to the 
quasiparticlc peak, one sees a broad maximum at larger 
Lo. This is a consequence of the IR lightcone. Over- 
ambitiously, if this were the correct interpretation, the 
location of the hump would give a measurement of the 
speed of light of the quantum critical theory. 
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APPENDIX A: SPINOR IN THE IR AdSi REGION 

The Dirac equation in the IR AdSii region including 
the mixing term is 

dr + L I Rinc^ / r 2icr^Lpori5LiR/r 



LiR f ka^ — ia'^uj/cm. 



^' = 



Now to solve this we employ the following basis rota- 



tion T2 
form: 



Then the Dirac equation takes the 



dr + —^{ka^ — ia'^uj/cm) 1 



L 



where 



IR 3 

— a ' 
r 



crljc-|(-fc,-w) 



* = (Al) 



(A2) 



We can now block diagonalize this equation into two in- 
dependent Dirac equations. We make the following basis 
rotation: 



-Vc 

+K 



Here, 



(A3) 



Vr = L 



IR 



u 



-2iijl^plLjii ~2ii%LplLjR 



(A4) 



where Uc determines the conformal dimension of the 
spinor in the IR AdS^ region. These Dirac equations are 
then exactly that of a spinor in AdSi with mass -izVc/ Lib. 
The (two) general incoming solutions can be found, and 
at the boundary of this IR AdS'4, a basis for these solu- 
tions behaves like 



Tl{k,uj) 
a^Tf{-k,-~cu) 



giR{k,uj)r--^ 





(A5) 



Tl\k,u;) 

a^Tf\-k,-u) 



( 







Qin{k,~uj)v-''^ 



(A6) 



V 



where the IR Green's function for a spinor is 



GiRik, 



r(l/2 - ly,) k + uj/cm f ,^2 
r(l/2 + yc) \ k- w/ciR 



(A7) 

We can then integrate these solutions out to the UV 
boundary where we can use similar methods to ([l3|) to 
read off a general form for the full Green's function. The 
result is (l46l). 
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